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U; ABSTRACT 

Qli Some recent observations seem to disagree with hierarchical theories of galaxy 

formation about the role played by major mergers in the late buildup of massive E-SO's. 
We re-address this question by analysing the morphology, structural distortion level, 
and star formation enhancement of a sample of massive galaxies (M* > 5 x 1O 1O M0) 
lying on the Red Sequence and its surroundings at 0.3 < z < 1.5. We have used an 
initial sample of ^1800 sources with K s < 20.5 mag over an area ~ 155 arcmin 2 on the 
Groth Strip, combining data from the Rainbow Extragalactic Database and the GOYA 
Survey. Red galaxy classes that can be directly associated to intermediate stages of 
major mergers and to their final products have been defined. We report observational 
evidence of the existence of a dominant evolutionary path among massive red galaxies 
at 0.6 < z < 1.5, consisting in the conversion of irregular disks into irregular spheroids, 
and of these ones into regular spheroids. This result implies: 1) the massive red regular 
galaxies at low redshifts derive from the irregular ones populating the Red Sequence 
and its neighbourhood at earlier epochs up to z ~ 1.5; 2) the progenitors of the 
bulk of present-day massive red regular galaxies have been disks that seem to have 
migrated to the Red Sequence mostly through major mergers at 0.6 < z < 1.2 (these 
mergers thus starting at z ~ 1.5); and 3) the formation of E-SO's that end up with 
M* > 10 n M© at z — through gas-rich major mergers has frozen since z ~ 0.6. All 
these facts support that major mergers have played a dominant role in the definitive 
buildup of present-day E-SO's with M* > 10 11 M Q at 0.6 < z < 1.2, in good agreement 
with hierarchical scenarios of galaxy formation. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: 
formation — galaxies: interactions — galaxies: luminosity function, mass function — 
galaxies: morphologies 
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1 INTRODUCTION 

Studies based on data from the Sloan Digital Sky Survey 
(SDSS) have revealed the existence of a co lor bimodality 
in the mass distribution of local galaxies (jStrateva et al.1 



* Affiliations are listed at the end of the paper, 
f E-mail: mpm@iac.es 



l200ll:lKauffmann et al.ll2003al : lBaTdrv et al.ll2004n . The most 
massive systems (basically spheroids) accumulate into a 
well-defined Red Sequence in color-magnitude diagrams, 
while the less massive blue ones (mos tly disks) reside into 
a more spread Blue Galaxy Cloud (j Strateva et all 1200 ll : 
iKaufrmann et al.ll2003bl ; iBaldrv et alj|2004 ). Although this 
color-mass bimodality is observed up to z ~ 1, observational 
data evidences a strong evolution in it for both field and 
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cluster environments ([Couch et al ]| 19981; Ivan Dokkum et al.l 



l200d ; lvan Dokkum fc Franxll200ll ; lBerl et al.ll2004bl) . In fact? 
the mass limit isolating the Red Sequence from the Blue 
Cloud rises as the redshift increases and the stellar mass 
harbored by the former has nearly doubled since z ~ 1 
(whereas that of the Blue Cloud has remained nearly con- 
stant , seelArnouts et al.ll2007l ; lLotz et alj|2008al : IStutz et"al] 
|200S| ; iTavlor et all l2009h . This points to a progressive 
buildup of the Red Sequence during the last ~ 9 Gyr, as- 
sociated with the migration of disks from the Blue Cloud 
to the Red Sequence through mechanisms that are still 
poorly understood, but that must be responsible of their 
star formation quenching and morphological tr ansforma- 
tion l|Brinchmann fc Ellisfeoocl ; iFaber et al.ll2007l , F07 here- 
after). 

According to hierarchical models of galaxy formation, 
the mechanism governing this evolution in the most mas- 
sive systems have been maj or mergers (i.e., with mass 
ratios from 1:1 to 1:4, s ee ISomerville fc Primackl 1 19991 : 
ISteinmetz fc Navarrd 120021 ). Present-day massive spheroids 
(E-SO's) are expected to be the result of the most mas- 
sive and violent merging sequences in the Universe, also 
being the latest systems to be completel y in place into 
the cosmic scenario (at z < 0.5, see |Pe Lucia et alj 
l200rj ; iHopkins et al.|[2008bl : lOser et alj|2010h . However, this 
prediction conflicts directly with recent data indicating 
that massive galaxies seem to have been in place be- 
fore their less-mass ive counterparts (a phenomenon known 
as downsizing, see Bundv et al. 2006; Cimatti e t alj [20061 : 



iPerez-Gonzalez et al.ll2008al ; iMortlock et alJkoill ). 

Also, this assembly epoch seems to depend 



the 



seems 
and its 



environment 



dThomas et al. 20051; di Serego Aliehieri et al. 20051. 


2006alfbl; IPerez-Gonzalez et al.l l2008bl; 


Cooper et al.l 


20101; iNiemi et al. 20ld; Vikram et al.l 2010 


). The wide 



range of ages found for E-SO's (~ 2-15 Gyr) and the 
disagreement between optical and NIR age estimators 
(sometimes, of up to ~ 6 Gyr for a given galaxy) indicate 
that E-SO's have been built up at different epochs and 
through differ ent mechanisms, b a sically depending on 
their masses dTrager et al.l |2000| ; iBregman et~ai1 |2006| ; 
lLopez-Corredoiral |2010| ~ It is generally accepted that the 
time period at 1 < z < 2 is a transition era in which an 
increasing fraction of galaxies end their star formation 



activ ity and move o nto the Red Sequence (lArnouts et al 



2001 but see also lEllid Jl997l; ICristobal-Hornillos et al 



2003; Elichc-Moral et al 



2006rJ 



ITavlor et al 



Brammer et al.l l2011al ; iDommguez Sanchez et al 



Nevertheless, there are conflicting views on the number 
evolution experienced by massive E-SO's at z < 1, with 
studies reporting from a negligible evolut ion to an increase 
by a factor of up t o ~ 3 in these systems llllbert et al.l 200C 



Faber et al 



Nicol et al. 



120071 ; IScarlata et all 120071 ; ICool et all 1200 
2011) 



Moreover, the Red Sequence is not made of a ho- 
mogeneous galaxy population, but of a mixing of dif- 
ferent ga laxy types t hat has evolved strongly with 

redshift l|Cimatti et al.l |2002| ; lYan fc Thompson! 120031 ; 



Gilbank et al.ll2003l; iMoustakas et al.ll2004l ; iFranzetti et all 
20071 ; iHempel et al .1 [20 1 lh . In the last years, different stud- 



ies have analysed the number evolution of the Red Sequence 
from color- and morphologically-selected samples, but the 



latter ones in general ( 


Franzetti et al.l 2007; Bundv et al.1 


2010i;lHempel et al.l2011 


). Recently. Ilbert et al.1 ( 2010T) have 



considered simultaneously both the morphological and star- 
formation properties of red galaxies to analyse their buildup. 
These authors confirm that the bulk of massive quiescent 
galaxies is rapidly created at 1 < z < 2, this mass assembly 
becoming negligible at later epochs. As in many previous 
studies, they also propose major mergers as drivers of this 
buildup. 

Several estimates assuming that most of the mass 
growth in quiescent galaxies is due to mergers indicate that 
this mechanism is capable of explaining at least 50% of 
the number density evolu tio n of massive galaxies (s ee, e.g., 
lEliche-Moral et al J lioiOal lbl; iBrammer et al.l l2011bl ). How- 
ever, no direct observational evidence has been found up to 
the date on the existence of a cause-and-effect link between 
major mergers and the definitive assembly of massive qui 
escent galaxies llFaber et al1l2007l; ISobral et al.1 [20091. 12011 



Bernard! et~aT1 l201Ct l201lT iLin et al.1 l20ld ; IChuter et al.1 



201 ll ) 



The present study tries to advance in the understand- 
ing of the formation and evolution of massive galaxies by 
analysing the physical properties of red galaxies at 0.3 < 
z < 1.5 with stellar masses M* > 5 x 1O 1O M . The nov- 
elty of our study over previous ones is two- fold. First, we 
include information about the structural distortion of each 
galaxy to trace merger remnants (besides considering mor- 
phology and star formation properties). And secondly, as 
most objects in their evolution towards the Red Sequence 
must have gone over nearby Green Valley locations transi- 
torily (F07), we have analysed the red galaxies both lying on 
the Red Sequence and at close positions on the Green Valley. 
Therefore, red galaxies in the context of this paper include 
both the galaxies on the Red Sequence and at its neighbour- 
hood. The galaxy classes resulting from the combination of 
morphological, structural, and star-formation activity prop- 
erties allow us to trace the evolution of intermediate stages 
of major mergers and of their final remnants since z ~ 1.5. 
Finally, the observed number density evolution experienced 
by each galaxy type is used to carry out a set of novel ob- 
servational tests defined on the basis of the expectations of 
hierarchical models, which provide observationally and for 
the first time main evolutionary paths among the different 
red galaxy types that have occurred in the last ~ 9 Gyr. 

The paper is organized as follows. In jj2l we provide a 
brief description of the survey. Section ij3] is devoted to the 
definition of the mass-limited red galaxy sample. In £01 we 
define the galaxy classes according to the global morphology, 
structural distortion level, and star formation enhancement 
of the red galaxies. In iJS] we comment on the sources of er- 
rors and uncertainties. Section fj6] presents three novel tests 
to check the existence of any evolutionary links between the 
different red galaxy types, based on the expectations of hi- 
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Figure 1. Color-color diagram for distinguishing stars from 
galaxies (/ — K s vs. V — I). Dots: Data from the original K- 
band selected catalogue. Solid line: Color cut defined to isolate 
stars (i.e., the sequence of data located below the line) from galax- 
ies (data lying above it). Circles: Objects classified as "stars" or 
"compact" visually (see All objects identified as "stars" 

visually are located in the stellar region of this diagram. 



Figure 2. Spectroscopic redshifts vs. photometric redshifts for 
the red galaxies in the sample having sp e ctrosc opic confirmation 
in the DEEP2 catalog l lDavis et alj|2003l . I2007T) . Bright and faint 
galaxies in each one of the three wide redshift bins under con- 
sideration in the study are plotted with different symbols. The 
typical photometric redshift uncertainty for the red galaxy sam- 
ple is below A(z)/(1 + z) < 0.03 at all redshifts for both bright 
and faint sources (see the estimates in the figure). 



erarchical models of galaxy formation. The results of the 
study are presented in SJ7] In particular, the results of the 
three tests proposed for the hierarchical scenario of E-SO for- 
mation can be found in i ]7.3l The discussion and the main 
conclusions of the study are finally exposed in §3H1 and |U 
respectively. Magnitudes are provided in the Vega system 
throughout the paper. We assume the concordance cosmol- 
ogy (O m = 0.3, Q.a = 0.7, and H = 70 km s" 1 Mpc" 1 , see 
ISpergel et aLll2007r i. 



2 THE SAMPLE 

We have combined multiwavele ngth data f r om th e Rain- 
bow Extragalactic Databas^B (| Barro et all 201 lal b |) and 
the GOYA photometric surve^f 2 ! (see iBalcells et al.l I2002T I 
over an area of 

Ah 



155 arcmin in the Groth Strip (a 



14 n 16 m 38.8 3 and S = 52° 16' 52", see iGroth et al.l Il994 
iRatnatunga et all 1 19991 ; ISimard et ail I2002T I. The~ltainbow 
Extragalactic Database compiles multi-wavelength photo- 
metric data from the UV to the FIR (and, in particu- 



1 Rainbow Extragalactic Database: 

https: / /rainbowx. fis.ucm.es/Rainbow_Database/Home. html 

2 GOYA project (Galaxy Origins and Young Assembly): 
http: / /www. astro. ufl.edu/GOYA/home. html 



lar, in Spitzer/MIPS 24[im band) over this sky area, pro- 
viding analysis of spectral energy distributions of nearly 
80,000 IRAC 3.6+4.5 /im selected galaxies. This study 
considers the photometric redshifts available in the Rain- 
bow Database, which have a typical photo metric red- 
shift a ccuracy of < Az/(1 + z) >= 0.03 (|Barro et all 
2011b), as derived for the sources with spectroscopic red- 
shifts available in the DEEP2 Galaxy Redshift Survey 
l|Davis et al.l 120031 . 120071 ). The GOYA Survey is a survey 
covering the Groth Strip compiling photometry in four 
optical bands (U, B, F606W, and F8UW) and in two 
near infrared ones (J and K s ) with visual classifications 
available, reaching similar depths to the Rainbow data in 
similar bands (U ~ 25 , B ~ 25.5, K ~ 21 mag, see 
Cristobal-Hornillos et al.l 120031: lEliche-Moral et all l2006bl: 



Abreu et al.ll2007i ; iDommguez-Palmero et alj|2008l ) 



We have performed the sample selection starting from 
a if-band selected catalog in the field, reaching a limiting 
magnitude for 50% detection efficiency at K ~ 21 mag. Sev- 
eral cuts have been performed to the original catalogue to 
obtain a mass-limited red galaxy sample. Firstly, red galax- 
ies are selected as detailed in §31 This selection determines 
the redshift interval of the study, as it is restricted to the 
redshifts where the obtained number of red galaxies is sta- 
tistically significant (i.e., to 0.3 < z < 1.5, see i|4.3|l . 
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According to the Mk-z distribution of the red galax- 
ies sample, the faintest absolute magnitude for which the 
catalogue is complete in luminosity up to z ~ 1.5 corre- 
sponds to Mjf,iim ~ —24 mag. According to the redshift 
evolution of the mass-to-light re l ation (assuming a Salpeter 
IMF) derived bv lArnouts et all |2007l ) for a sample of qui- 
escent bright galaxies, a red passive galaxy with this K- 
band absolute magnitude at z ~ 1.5 has a stellar mass of 
A/*,iim ~ 5 x 10 10 M Q . Therefore, we have selected red galax- 
ies with masses higher than M* ~ 5 x 1O 1O M0 at each z just 
considering all galaxies with Mjc,cut(*0 = —23.3 — 0.452 to 
account for their luminosity evolutio n. This luminosity cut i s 
very similar to the one obtained bv lCirasuolo et all |2007l y 
Analogously, the equivalent luminosity cut to obtain a com- 
plete red galaxy sample for M*,u m = 1O 11 M0 would be: 
Mk = —24 — 0.452 (we will use this selection for compar- 
ing our results with those reported in other studies, see JJ}. 
The detection e fficiency in the Jf-band drops b elow 0.9 for 
ma > 20.5 mag (ICristobal-Hornillos et al1l2003l ). so we have 
checked that all galaxies in our mass-limited red sample ex- 
hibit apparent magnitudes brighter than this limit. 

We have used the color-color diagram shown in Fig.[T] 
to remove stars from the sample. It represents the I — K s 
vs. V — I distribution for all the sources in the mass-limited 
red galaxy sample. Stars typically exhibit NIR colors bluer 
than galaxies, so they populate the lower region in the di- 
agram. Attending to this bimodality of star-galaxy colors, 
we have defined a color-color cut to isolate galaxies from 
stars (see the solid line in the figure). The marked points in- 
clude the stars and compact objects in the sample. We have 
checked that the stars identified according to it include all 
the objects at lower redshifts that have been classified as 
"stars" in the morphological classification performed in tj4.ll 
(they are marked in Fig.[T]). As the number of compact ob- 
jects found in the sample is not statistically significant (see 
H4.ll) . they have been excluded from this study. From an 
initial sample of 1809 sources from the original Jf-band se- 
lected catalogue, we finally have a mass-limited red galaxy 
sample of 257 systems at 0.3 < z < 1.5. 

Figuxe[2] compares the photometric and the spectro- 
scopic redshifts for the red galaxies in our mass-limited 
sample with sp ectroscopic redshif t determinations from the 
DEEP 2 survey (|Davis et al.ll2007T ). The average redshift un- 
certainties are similar for both bright and faint sources at 
all redshifts, being below A(z)/(1 + z) ~ 0.03 in the whole 
redshift interval under study. 



3 RED GALAXIES SELECTION 

Previous studies on red galaxies have been traditionally 
centered on the galaxies lying on top o f the Red Se- 
quen c e (see, e.g., van Dokkum et al. 1 12000| ; IWillmer et al.l 
120061 ; ICirasuolo et al.l 2007 ) . These studies derive a linear 
fit to the rest-frame (U — -B)rcst-frame-Afs relation and con- 
sider all the objects within la scatter of this relation as 
red galaxies. Nevertheless, the majority of the galaxies that 
have migrated to the Red Sequence is expected to have 



spent part of their lives at nearby locations on the Green 
Valley, undergoing transitory stages of ~ 1 Gyr in their 
evolution towards typical E-SO's (F07). Many observational 
and computational studies confirm that advanced stages 
of major mergers usually lay on the Red Sequence or at 
nearby locations on the Green Valley, because gas-rich en- 
counters experience noticeable dust exti nction due to the 
starbursts induced by the encounter (seelKeel fc Wulll995l; 
Bell fc de Jon3|200l|; [kauffmann et al.ll2003bl; iRoche et al " 



20061; |Cox et al.ll2008l; ICirasuolo et al 

2009al ; IChilingarian et al.ll20ld ; lLotz et al.ll2008bl . l2010al lbh 



2007 



Poggianti et al 



Therefore, in this study we have included the systems both 
lying on the Red Sequence and on nearby positions on the 
Green Valley into our red galaxy sample, to also trace these 
transitory evolutionar y stages towards the Red Sequence 
(see more references in lEliche-Moral et al.ll2010af ) . 



3.1 Color cuts for red galaxies selection 

Following previous studies, we have selected red galaxies 
according to thei r rest-frame U — B color (see, e.g., F07; 
iKriek et all 120081 ). This color traces the 4000 A break up, 
a spectral feature characteristic of evolved stellar popula- 
tions. However, the rest-frame U — B estimates derived from 
the present dataset exhibit so high uncertainties due to 
photometry and redshift-determination errors that we have 
used instead the observed color that samples the rest-frame 
U — B more closely at the center of the three wide redshift 
bins covering the complete redshift range (0.3 < z < 0.7, 
0.7 < z < 1.1, and 1.1 < 2 < 1.5). 

The top panels of Fig.[3] show the color distribution of 
all galaxies in our Tf-band selected catalogue in the colors 
that sample more closely the rest-frame U — B at these three 
redshift bins (U — V , B — I, and V — K for the low, mid- 
dle, and high redshift intervals, respectively). Typical color 
errors of the galaxy sample are indicated with a horizontal 
bar at the left top corner of each frame. The figure shows 
the well-known bimodal distribution of galaxies into red and 
blue populations at all redshifts up to 2 ~ 2 (see references 
in EJXJ) . These bimodal distributions have been modelled as 
the addition of two Gaussian functions (also plotted in the 
figure) . 

In order to include the galaxies located at nearby posi- 
tions on the Red Sequence in the red sample (and not just 
the galaxies lying within la scatter of the red peak), we have 
adopted the colors at which both Gaussian distributions 
cross as the color cuts for isolating red from blue galaxies at 
each redshift bin. The resulting cuts are (U — V) = 1.4 mag 
for 0.3 < 2 < 0.7, (B - I) = 2.7 mag for 0.7 < 2 < 1.1, and 
(V-K) = 4.5 mag for 1.1 < 2 < 1.5. 

In the bottom panels of Fig.[3l we show the histograms 
of the rest-frame U — B color in the same redshift bins. 
The dashed red and dotted blue histograms show the dis- 
tribution of red and blue galaxies selected according to the 
criteria based on the apparent colors commented above. In 
1.1 < 2 < 1.5 redshift range, the red and blue galaxy distri- 
butions appear quite mixed, probably due to the high errors 
associated to rest-frame U — B colors (see the horizontal bar 
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Figure 3. Color histograms of all galaxies in the X-band selected catalogue in the three wide redshift bins under consideration (0.3 < 
z < 0.7, 0.7 < z < 1.1, and 1.1 < z < 1.5). Top panels: Histograms of the apparent color nearest to the rest-frame U — B at each redshift 
interval. Red and blue solid lines: Gaussian functions fitted to match the bimodal color distributions at each redshift. Vertical dashed 
lines in top panels: Color cuts used to isolate red from blue galaxies in each redshift bin (see the text). Bottom panels: Histograms of 
the rest-frame U — B in the same redshift bins. The horizontal bar in the upper left corner of each frame represents the average color 
error of one galaxy in the sample in each case. Dashed red and dotted blue histograms in these panels show the distributions of red and 
blue galaxies, respectively, selected according to the criteria based on the apparent colors of the top panels. Red galaxies are selected on 
the basis of apparent colors just attending to the low uncertainties associated to them compared to the typical errors of the rest- frame 
U — B colors at all rcdshifts. 



in the upper left corner of each frame). These high errors in 
the rest-frame U — B color inhibited us to perform the red 
galaxies selection on the basis of rest-frame U — B colors. 

Previous studies report a negligible dependence of the 
rest-frame U — B color cut isolating the Red Sequence from 
the Blue Cloud with the galaxy mass for both the mass and 
redshift range considered in t his study (the U — B color cut 
varie s ~ 0.3 mag at most, see I Taylor et al. 2009; Nico l et al.l 
2011). Therefore, any dependence of our color cuts with mass 
(or, equivalently, with luminosity) has been overridden. 

3.2 Galaxy populations in the red galaxy sample 

We have analysed the kind of galaxy populations selected 
at each redshift interval with the color cuts defined in £|3 . 1 1 
to find out if the selection is homogeneous at each redshift 
bin. In Fig. U we show the color-redshift distributions of all 
the galaxies in the if-band selected catalogue, for the color 
indices nearest to the rest-frame U — B color in the three 
wide redshift bins considered in the study (0.3 < z < 0.7, 



0.7 < z < 1.1, and 1.1 < z < 1.5). The color cuts to distin- 
guish between red and blue galaxies at each redshift bin are 
marked in their corresponding panels with horizontal lines. 
The red galaxies selected at each redshift bin are those above 
the color cut and enclosed between the redshift limits asso- 
ciated to the color index in each panel. 

We have overplotted the theoretical evolution followed 
in each color-redshift plane by different galaxy types. 
These tren ds have been model led using the IRAF package 
C0SM0PACK (|Balcells et al.ll2003l ). which uses the spectral en- 
ergy distributions (SEDs ) predicted by the s t ellar popula- 
tion synthesis models by iBruzual fc Charlotl (120031 ) to ob- 
tain the evolution of color indices with z. The galaxy types 
plotted in the figure include E's and spirals formed at dif- 
ferent redshifts, a SO, a star-forming irregular, and a dust- 
reddened star-forming galaxy. Standard star-formation his- 
tories (SFHs), metallicities, and dust extinction values are 
assigned to each galaxy type, accordingly to observations 
llKennicutdll989l ; iGallazzi et "aLll2005l ; iMufioz -Mateos et alj 
2009). The modelling parameters assumed for each simu- 
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Table 1. Simulated star formation histories characteristic of different galaxy types. 



Galaxy Type 


SFR law 


Metallicity 


r»,v (mag^ 


Zf 


(1) 


(2) 


(3) 


(4) 


(5) 


Elliptical 


Finite burst (r = 0.7 Gyr) 


0.02 


0.6 


0.8,1.0,1.3,1.5,2.0,2.5,3.0 


SO and spirals 


Exponentially declining with 
e f = 1, 2, 3, 5, and 7Gyr^ 


0.02 


0.6 


1.2,1.3,1.5,2.0,3.0 


Irregular 


Constant 


0.008 


0.0 


Present at all z 


Dust- reddened 
star-forming galaxy 


Constant 


0.008 


3.0 


Present at all z 



r ti ir: Optical depth of the stellar population in the V-band due to dust extinction. 
ef. Star formation e- folding timescales in exponentially declining SFRs. 




0.5 1 1.6 2 2.5 3 0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

z z z 

Figure 4. Color-redshift distributions for all galaxies in the i\-band selected catalogue, in the observed colors closest to the rest-frame 
U — B color in the three wide redshift bins of the study. Panel a: U — V color (used for red galaxy selection at 0.3 < z < 0.7). Panel b: 
B — I color (used for red galaxy selection at 0.7 < z < 1.1). Panel c: V — K color (used for red galaxy selection at 1.1 < z < 1.5). Dots: 
Observational data. Lines: Theoretical trends in the color-2 space for different galaxy types according to standard SFHs and physical 
properties of different galaxy types (more details in the text). Red: Ellipticals with different Zf . Magenta: SO galaxy with Zf = 3. Yellow: 
Spirals with different Zf. Blue: Star-forming irregular. Green: Dust-reddened star-forming galaxy with r* y = 3. Horizontal black lines: 
Color cuts defined to isolate red from blue galaxies, according to Fig. [3] Vertical solid lines: Limits of the redshift interval associated to 
the color plotted in each panel. Red galaxies in each redshift bin are those enclosed between the vertical lines and above the horizontal 
line marking the color cut. 



lated galaxy type are listed in Table[T] We must remark 
that we have adopted for the ellipticals a finite burst in- 
stead of the traditional single stellar population model. Al- 
though both SFRs provide similar colors and magnitudes 
evolution, the finite burst model can account for the short 
initial phases of st ar formation observed in these sys tems 
(|Kannappan et a l. 2009; Hucr tas-Companv et al.ll2010l ). 

Different formation redshifts have been considered for 
all the types. Nevertheless, as the redshift evolution of colors 
for SO's and their values were quite independent of Zf, we 
have just plotted the SO track for zt — 3 in Fig.|4]in benefit of 
clarity. As Irr's with different zt are basically located below 
the spirals region in the plots, we only plot the case with 
Zf = 3 for the same reason. We have also simulated several 
cases of dust-reddened star-forming galaxies with different 
dust extinction levels. Their tracks basically overlap with 



those of the ellipticals in the plots, so we only show the 
model with T*,y = 3 mag. 

Although we have assumed different Zf for the forma- 
tion of spirals, to assume the same epoch for their formation 
have a negligible effect on both their color evolution and lo- 
cation in the color- z plane. They differ basically in their e- 
folding timescales, metallicities, and gas and dust contents. 
However, a similar conclusion cannot be derived for ellip- 
ticals. If we assume that all ellipticals have formed at an 
early epoch (e.g., Zf ^ 2, as supported by some studies, see 
Jl]), the observed color- z space of real red galaxies cannot 
be reproduced with the tracks modelled for Zt > 2.0, unless 
there is a large population of dust-r eddened star-forming 
objects at all redshifts at 0.4 < z < 2 faitzbichl er fc White] 
120061 ; lEliche- Moral et aI1l2010al ). A wide range of values for 
the formation redshift of ellipticals reproduces more prop- 
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erly the observed color-z distribution of red galaxies (see 
the mesh of red lines in each panel of the figure), a fact 
that does not exclude the existence of relevant populations 
of dust-reddened objects at different redshifts. 

The models for E's, SO's, and dust-reddened star- 
forming galaxies lie all above the color cuts in the three 
panels and that no spiral track enters the region of red 
galaxies in Fig.U This proves that the color cuts used in 
this study select quite homogeneous samples of galaxy pop- 
ulations in the whole redshift interval, basically galaxy types 
earlier than Sa and dust-reddened star-forming galaxies. 

3.3 Comparison with other studies 

The red galaxy selection made in this study cannot be di- 
rectly compared to the red galaxy samples obtained by most 
studies in the literature because, first, we have included red 
galaxies adjacent to the Red Sequence to study objects at 
transitory stages of their evolution towards it (which is not 
usual, see JTJ, and secondly, we ha ve estimated m a sses u sing 
the M 4 /Lk-z relation derived by lArnouts et al.l l|2007h for 
different redshifts, wher eas most autho r s use the M,/Lv- 
color relation derived bv lBell fe de Jonel l|200lf ) or an equiv- 
alent relation. Moreover, most studies report the number 
evolution of red galaxies for masses M* > 10 1 Mq, instead 
of for masses M* > 5 x 10 10 Mq (as our case). In order to 
check out our results, we have made alternative red galaxy 
selections for M* > 10 n Mo, adopting the color cuts and/or 
the mass-to-light relation used by other authors. 

The three panels of Fig. [5] compare the redshift evolu- 
tion of the number density of red galaxies derived from our 
data with the results of different authors, for analogous mass 
and color selections in each case. In pa nel a, we have assumed 
the U — B color evolution derived by Ivan Dokkum fe Franxl 
(2001) to select red galaxies (following F07), and the 
masses are est imated using the M./Lv-color relation by 
iBell fc de Jond . Only red galaxies with M* > 1O 11 M at 
each redshift are considered in this panel. Panel b of the 
figure also assumes the color cut by Ivan Dokkum fe Franxl 
for selecting red galaxi es, but the ma ss estimates assume 
the M*/L/f relation by lArnouts et all , which includes evo- 
lutive corrections (it is equivalent to the one derived by 
ICirasuolo et~ai1l2007l ). The number densities of red galaxies 
shown in panel b are also for galaxies with M, > 10 1:L M© 
at each redshift. Finally, panel c of the figure use the same 
selection criteria as panel b, but the number densities of red 
galaxies have been computed for M„ > 5 x 1O 1O M0. Note 
that the results of our color selection (including galaxies on 
the Red Sequence and at nearby locations) are not plotted 
in this figure (see $7] and Tabled . 

The data frorn lArnouts et all (|2007r ). ICirasuolo et al.1 
120071 ). and F07 have been obtained by integrating their 
red galaxy luminosity functions at each redshift for 
log(A1»/MQ) > 11. In this case, the absolute magnitudes 
have been transformed into stellar masses using t he ex - 
pression derived for local E-SO's by ICimatti et all (2006), 
considering the L-evolution of red galaxies derived by F07 
and the redshift evolution of the (B — K) color expected 



for E - SO's JShimasaku fc Fukugital 19981 ; Ide Lapparent et al.l 
120031 ; lArnouts et al-lbOolT AB magnitudes have been trans- 
formed to_the_yeg_a_sjr^mJrithe B and K bands accord- 
ing to iBlanton fc Roweisl (|2007f ) transformations and con- 
sidering galaxy colors derived for E-SO's by iFukugita et all 
|l995l ). 

The good agreement between our results and those from 
independent studies for similar selection criteria supports 
the reliability of our methodology and completeness of our 
nominal red sample (compare the black filled circles with the 
rest of studies in all panels of Fig.(5]). However, although our 
first data point in panels b and c in the figure is inside the 
cloud of points within errors, it does not follow the trend of 
the other authors. So, our data at this redshift are probably 
affected by volume and cosmic variance effects. 



4 CLASSIFICATION OF RED GALAXY 
TYPES 

The hierarchical picture of galaxy formation predicts that 
massive E-SO's are the result of the most violent and mas- 
sive merging histories in the Universe (see references in §1 in 
EM10). To test this scenario, we need to distinguish between 
galaxies undergoing a major merger and normal E-SO's (see 
jJBJ. Normal relaxed galaxies and major mergers differ basi- 
cally in their structural distortion level. Major mergers also 
exhibit different global morphology and star formation en- 
hancement depending on the gas content of the progenitors 
and the evolutionary stage of the encounter. 

A gas-rich major merger is expected to turn into a dust- 
reddened star-forming disk with noticeable structural distor- 
tions at intermediate and advanced stages of the encounter, 
basically since the coalesce of the two galaxies into an unique 
galaxy body (this stage is known as the merging-nuclei 
phase). In earlier phases of the merger, the two galaxies 
can develop noticeable tidal tails and asymmetric structures, 
but the two bodies can still be distinguished and are not ex- 
pected to suffer from enough dust reddening to lie nearby 
or on top of the Red Sequence. During the latest phase of 
the encounter (post-merger), the star formation is quenched 
and the remnant gets a more relaxed sph eroidal structure 
until it transform s into a typical E-S0 fsee lCox et al.ll2008l ; 
lLotz et al.1 l2008bl . l2010al lbf). Intermediate-to-late stages of 
gas-poor major mergers present a distorted spheroidal mor- 
phology and n egligible levels of st ar formation, thus being 
quite red too (|van Dokkumll2005h . By the contrary, typi- 
cal E-SO's present a spheroidal-dominated relaxed morphol- 
ogy, although they are also expected to be quite red due to 
their negligible star formation. Therefore, a gas-rich major 
merger is expected to be quite red from its merging-nuclei 
phase until its transformation into a typical E-S0, so we have 
traced major mergers once the two merging galaxies have 
merged into a unique remnant, because they are expected 
to be quite red in any case. Accounting for this, we have 
classified the galaxies in our red sample attending to their 
global morphology, structural distortion level, and star for- 
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Figure 5. Comparison of the redshift evolution of the number density of red galaxies obtained in this study with results from the 
literature for different mass and color selection criteria (see the text for details). Panels a and b plot the results for red galaxies with 
M» > lO n M , but for different mass-to-light transformation relations. Panel c presents the results obtained with the same selection 
criteria as panel b, but for red galaxies with M* > 5 X 10 10 Mq. All data plotted in the same panel use equivalent color and mass selection 
criteri a . Filled black circles: R esults for our if-band sel e cted data for each s election. Rest of symbols: R esults obtained bylArnouts et alj 
j2007l ), lc7rasuolo et af]<2007ft . F07. lTavlor et al.l fl2009t) , IHbert et alj fcOld l. lMatsuoka fe Kawaral J2010T ), and lNicol et all fcoilll . Plotted 
data assume h = = 0.7 and a Salpeter initial mass function. Cosmic variance and sample incompleteness contribute to the large 

dispersion found among different studies at z < 0.4. 
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Figure 6. False-color postage stamps of some red galaxies in our sample, obtained using the V and I bands. One example representative 
of each type is shown for each wide redshift bin used in the present study (0.3 < z < 0.7, 0.7 < z < 1.1, and 1.1 < z < 1.5). North is up, 
East is left. The frames correspond to a 5" X 5" field-of-view, except for the irregular spheroid at 0.3 < z < 0.7 (141618.91+521352.3), 
where a 10" X 10" view is used to emphasize the interacting group to which this galaxy belongs (the galaxies at the North and at the 
Southern-East of the central object are located at the same redshift as the central object). No compact or unclassified objects are found 
at 0.3 < z < 0.7. 
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Figure 7. /-band contour maps and surface plots of some red galaxies in Fig. [6] One galaxy representative of each type is shown ("C", 
"RS", "RD", "IS", "ID", "NC"). The axes in both the contour and surface maps are in CCD pixels. 



mation enhancement to distinguish among normal galaxies 
and intermediate-to-advanced phases of major mergers. 

4.1 Classification according to global morphology 
and distortion level 

Despite the de finition of very effic ient quantitative morpho- 
logical indices (|Lotz et al.ll2008al ). visual inspection is still 
the m ost trustworthy method to classify galaxies morpholog- 
ically (|Schawinski et al.ll2007l ; IJogee et al.ll2009l ; iDarg et al.1 
120101 ). The red galaxies in our sample have been thus clas- 
sified visually by three co-authors. The comparison of the 
visual results with those obtained through quantitative clas- 
sification methods has proven the reliability and robustness 
of the visual classification (see i)4.2[) . 

No preconceived or classical morphological types have 
been used for the visual classification, as the emer- 
gence of the classi cal Hubble typ e s seems to occur at 
z ~ 1.0-1.5 dOdewahn et all 1 19961 ; iPapovich etafl 120051 ; 
iRavindranath et al.l 120061 ). The visual types of red galax- 
ies have been defined attending to the global characteristics 
exhibited by these galaxies in the sample. We find that red 
galaxies in our sample can be grouped morphologically at- 
tending to two aspects: 1) their structural distortion, and 2) 
its disk- or spheroid-dominated morphology. Six major ex- 
clusive classes have been identified according to the previous 
two criteria (see Fig.[6|: 

(i) Compact galaxies (C).- Galaxies exhibiting compact 
morphologies, according to the seeing of the /-band images. 
No spatial information is available for them. The number 
of these objects in the sample is negligible, so we will not 
consider them henceforth. 

(ii) Regular Spheroids (RS).- Galaxies with regular 
isophotes and dominated by a central spheroidal component. 
We remark that our RS's do not correspond t o the spheroidal 
type defined bv lKormendv fc Bender! ([2012) . This class ba- 
sically groups E-SO galaxies. 

(iii) Regular Disks (RD).- Galaxies with regular isophotes 
and dominated by a disk component. Dust-reddened, typical 
spiral Hubble types correspond to this class. 

(iv) Irregular Spheroids (IS).- Galaxies with irregular 



isophotes in the whole galaxy body but dominated by an 
spheroidal component. 

(v) Irregular Disks (ID).- Galaxies with irregular 
isophotes in the whole galaxy body and dominated by a 
disk. 

(vi) Non-Classified or unclassified objects (NC).- Galax- 
ies that cannot be classified into any one of the previous 
classes, because of their faintness or noise. The number of 
these objects is also negligible, so we will not consider them 
from now on. 

We have used surface brightness isophotes and surface 
maps to identify structural or morphological features that 
are more noticeable in these maps than in normal images 
(see Fig. [7]). We have also adopted the following additional 
classification rules. As commented above, we have traced ex- 
clusively advanced stages of major mergers, i.e., once both 
galaxies have merged into a unique body. Therefore, we have 
assigned a certain type to a red galaxy regardless of its en- 
vironment, i.e., independently of whether it has close neigh- 
bours or not. We have considered as irregular galaxies only 
those systems that exhibit a noticeable distortion level in 
its whole body. This implies that a galaxy in an interacting 
pair is identified as an independent regular galaxy, because 
its central body does not exhibit noticeable distortions yet, 
despite of the existence of significant tidal features in its out- 
skirts. We remark that we are excluding from the irregular 
class the early phases of mergers in which the interacting 
galaxies have not merged into one body yet. The reason is 
because the tests for the hierarchical scenario of E-SO forma- 
tion described in Sj6]are based in the tracing of these specific 
advanced stages of major mergers (and not of earlier phases). 

Minor mer gers imprint less significan t distortions than 
major mergers (jEliche-Moral et al.|[2006al . l201ll ; | Lotz et al.l 
l2008bl ). but they can also redden the galaxy enough to locate 
it a t neighbouring re gions of the Red Sequence temporar- 
ily l|Lotz et al.ll2010aH bl). We have avoided the inclusion of 
minor mergers into the irregular class by considering exclu- 
sively morphological features typical o f major mergers t o 
classify a system as irregular, following Ijogee et al.l |2009h . 
such as the existence of multiple nuclei of similar luminos- 
ity in the body, the existence of equal-length tidal tails, or 
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"train-wreck" morphologies. So, our red galaxy sample con- 
tain some minor mergers, but they are included into the 
regular class. 

Our classification criteria try to minimize the effects of 
the spatial resolution loss and cosmological dimming inher- 
ent to the rise of redshift, which soften the luminosity distri- 
bution of galaxies, losing the physical regions of fainter sur- 
face brightness. Therefore, we have defined our galaxy types 
attending to features that affect to the whole galaxy body 
and that are poorly affected by these effects for the consid- 
ered redshifts. Moreover, the classification has been done in 
the reddest visual band available with the best spatial reso- 
lution (/-band from HST/WFPC2). This band samples the 
rest-frame visual spectrum in the whole redshift range under 
study (from I at low redshifts to B at z ~ 1.5). Galaxies 
usually exhibit more distorted and clumpy morphologies in 
the rest-frame UV than in optical or NIR bands, but the 
morphology of a system is very similar in all rest-fra me op- 
tical bands |Kuchinski et al.ll2000l ; [jogee et al.ll2009t) . so we 
expect our visual classifications to not be biased towards ir- 
regular types at any redshift in our study. Simulations of the 
cosmological effects in the apparent /-band morphology of 
our systems have demonstrated the robustness of our visual 
classification against this effect (see i)4.2.2[) . In general, the 
three independent classifiers agreed in ~ 90% of the classi- 
fications at each redshift bin, a fact that strongly supports 
the robustness of the visual classification. 

4.2 Tests to the robustness of the visual 
morphological classification 

We have carried out two kind of quantitative tests to check 
the robustness of the visual classifications: 

(i) Tests on the reliability of the visual distinction 
of spheroid/disk-dominated and irregular/regular classes, 
based on the quantification of the light concentrations and 
the deviations of the isophotes from perfect ellipses in the 
whole galaxy body. 

(ii) Tests on the cosmological effects on the visual mor- 
phology in the observed /-band (resolution loss, cosmolog- 
ical dimming, and change of the observed rest-frame band 
with redshift). 

4-2. 1 Reliability of the visual classifications 

In order to check our visual classification, we have com- 
puted concentra tion and asymmetry indices as defined by 
IConselicd (|2003l 1. which are quantitative indices closely re- 
lated to our visual classes by definition (irregular galaxies 
should exhibit high asymmetries, while spheroid-dominated 
galaxies must have high concentrations). We find that galax- 
ies visually classified as RS's show high concentration indices 
in general (~ 90%). The irregulars (IS's and ID's) present a 
wide spread of concentration values, consistently with their 
merger-related nature. 

Nevertheless, we also find that the galaxies classified 
visually as irregulars do not exhibit any correlations with 



the asymmetry index, contrary to expectations. This is 
probably because the asymmetry index estimates are ex- 
tremely affected by background substructures, thus requir- 
ing high signal-to-no i se data to be reliable (S/N > 100, se e 
IConselice et~al]|2000l ; IConsehcdl2003l ; IConselice et~al1l2009l ) . 
But the galaxies in our red sample have S/N ~ 40 — 50 
at most (as /-band magnitude errors are ~ 0.02 — 0.03 mag 
typically), making our asymmetry estimates quite uncertain. 
Moreover, asymmetry indices are sensitive to environmental 
influences in the galaxy outskirts. This means that some 
galaxies identified as regulars according to our criteria (be- 
cause they do not exhibit noticeable distortions in its whole 
body) may have a high asymmetry index because of tidal 
features in the outer parts. This obviously smudges the cor- 
relation between visual irregularity and computed asymme- 
tries. 

W e have adapted the method by IZepf fc Whitmorel 
(|l993h to quantify the irregularity level of galaxy morphol- 
ogy. These authors developed a procedure to classify an el- 
liptical galaxy as regular or irregular, attending to the dis- 
tortion level of their isophotes with respect to perfect el- 
lipses. They fitted ellipses to the isophotes of each elliptical 
to obtain the radial profiles of the coefficients 03, 63, 04, 
and 64 of their Fourier expansion series. The peak value of 
each Fourier coefficient was identified along the radial pro- 
file. These authors considered the following criteria to dis- 
tinguish among regular and irregular galaxies: 

(i) If all the peaks of the coefficients were small, this 
meant that the isophotes exhibited small deviations from 
perfect ellipses. Therefore, these galaxies could be consid- 
ered as regulars. 

(ii) If one of them was not small, then they differentiated 
between two possible cases: 

• If the maximum value of the peaks of all coefficients 
did not correspond to 64, then the isophotes deviated no- 
ticeably from ellipses, meaning that the galaxy was irreg- 
ular. 

• If the peak of 64 was the maximum among all the 
peaks, its classification depended on its trend with the 
radial position in the galaxy. If the profile of this coeffi- 
cient changed from one type to another within 1.5 effective 
radius of the galaxy, then the galaxy was irregular. If not, 
it just implied that the galaxy was boxy or disky (depend- 
ing on the sign of 64), but the galaxy morphology could 
be considered as regular. 

We can adopt this method for our galaxies, as we are 
considering irregularities that must affect to the galaxy as 
a whole. We have limited the analysis in each red galaxy to 
the isophotes with a mean signal higher than 1.5 times the 
standard deviation of the sky signal per pixel. We have used 
the IRAF task ellipse for fitting ellipses to the isophotes 
and for getting the third and fourth order coefficients of 
their Fourier expansion series. We have identified the peaks 
of each coefficient in the galaxy radial profile. 

The task ellipse uses a normalization to the surface 
brightness of the isophote that directly measure the devia- 
tions of the isophote from perfect ellipticity. According to 
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Figure 8. Comparison of visual and quantitative morphological classifications into regular /irregular classes. The histograms of the 
galaxies corresponding to a given visual class are shown as a function of the quantitative irregularity parameter Cj rr) i SO ph defined in 
£14.2.11 for the three wide redshift bins of the study. Left panels: For 0.3 < z < 0.7. Middle panels: For 0.7 < z < 1.1. Right panels: 
For 1.1 < z < 1.5. Visual regulars, irregulars, and stars or unclassified objects are plotted from the top to the bottom panels of 
each column, respectively. The galaxies classified as regular/irregular according to the quantitative method are indicated in each panel. 
Vertical dashed lines: Limiting values of Cj rrj i sop h to classify a galaxy as regular (|Ci rr ,isoph| = 0.05). Galaxies with |Ci rr ,isoph| < 0.05 
are regular. However, | Cirr ,isoph I values above this limit do not necessarily imply irregular morphology (see the cases in the text). Visual 
and quantitative classifications of irregularity agree in ~ 75% of the cases at all redshifts on average, supporting the robustness of the 
visual method. 



de lPeletier etal] <jl99(f ). these deviations can be considered 
negligible if they are < 5%, a value that can be translated 
directly to a value of 0.05 in these coefficients. Therefore, 
we have adopted this limit as the critical value to distin- 
guish between small and high values. We have defined the 
irregularity index Ci rr ,isoph as the peak value of maximum 
absolute value among the peaks of the four Fouri e r coe f- 
ficients. Therefore, according to IZepf fc Whitmord (|l993h . 
any galaxy that has |Ci r r,iso P h| < 0.05 is regular. If not, it is 
irregular, except if Cj r r,iso P h corresponds to the &4 coefficient 
and this coefficient does not change between [&4| > 0.5 and 
1 64 1 < 0.5 values or viceversa along its radial profile. 

We compare the results of the visual and quantitative 
classifications concerning the irregularity level of our red 
galaxies in Fig.[8] for each wide redshift bin. The percent- 
age of agreement between the visual and quantitative clas- 
sifications into the regular type is ~ 77% (decreasing from 
83% to 70% from low to high redshifts). This percentage is 
slightly lower in the irregular type: ~66% (rising from 58% 
to 74% from low to high redshifts). The miss-classifications 
between both methods are ~23% for visually-identified reg- 
ular galaxies (rising the confussion percentage from 17% to 
30% from low to high z), and ~34% for visual irregulars 
(dropping from 42% at low z to 26% at high z). In general, 
both procedures coincide in ~ 78% of the galaxy classifica- 
tions at 0.3 < z < 0.7, in ~ 68% of the classifications at 
0.7 < z < 1.1, and in ~ 73% at 1.1 < 2 < 1.5. 

It is important to remark that only isophotal data above 



1.5 times the sky standard deviation have been considered 
to quantify the irregularity of the galaxy in the quantita- 
tive method, i.e., it analyses isophotes of S/N > 1.5 at all 
redshifts. This implies that the quantitative method is not 
biased towards more irregular types at high redshifts, as it 
is limited to the isophotes with enough S/N at all redshifts. 
Obviously, Ci rr ,j sop h is derived from an intrinsic physical re- 
gion in the galaxies smaller at high redshift than at low 
redshift, just because of cosmological effects. But, as com- 
mented above, we consider as irregular galaxies only the 
stages of advanced major mergers, which imply a noticeable 
distortion level in its whole body. The effects of cosmolog- 
ical dimming and resolution loss on the classification are 
analysed in 34.2.21 

In conclusion, this test proves the robustness of the vi- 
sual classifications into regular and irregular types at all 
redshifts. 



4-2.2 Robustness of the observed morphology against 
cosmological effects 

In order to find out how the loss of spatial resolution, the 
cosmological dimming, and the change of rest-frame band 
with redshift are affecting to our classification, we have simu- 
lated images of galaxies at different redshifts in the obse rved 
/-band. We have used C0SM0PACK (|Balcells et all 120031 '). an 
IRAF package that transforms images of real galaxies to de- 
pict their appearance at a given redshift as observed with a 
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given telescope, camera, and filter. The transformation in- 
cludes K-corrections, change of observing band, repixelation 
to the scale of the observing system, convolution by the see- 
ing, and noise from sky, detector, and dark current. 

Starting from the /-band image of a galaxy represen- 
tative of one type at the lowest wide redshift bin (0.3 < 
z < 0.7), we have simulated the observed /-band image of 
the same object at the middle of the other two redshift bins 
(z — 0.9 and z — 1.3), in order to check if its morpholog- 
ical classification changes as the object is placed at higher 
redshifts. We have assumed a typical SFH, metallicity, and 
formation redshift for the galaxy to assign a characteristic 
SED at its original redshift. Once we obtain the simulated 
image at z = 0.9 and z — 1.3, we have classified them visu- 
ally following the same procedure as with the original image. 

The surface maps at different redshifts show that the 
distinction between the spheroid/disk-dominated morphol- 
ogy is quite robust in the whole redshift interval under con- 
sideration. In general, we can distinguish the winged surfaces 
of disks in the surface maps at all redshifts (see Figs.|6][7]). 
This means that the visual classification into spheroid/disk- 
dominated systems is quite robust against cosmological ef- 
fects for the redshifts and magnitudes studied here. 

Concerning the irregularity, the contour plots support 
the expected trend of external isophotes to appear more ir- 
regular as we move towards higher redshifts. But this ef- 
fect should not affect strongly our results, as we derive 
our visual classification into regular/irregular classes hav- 
ing into account the morphology of the isophotes over the 
whole galaxy body, discarding the most external ones. This 
qualitative test indicates that cosmological effects are not 
expected to affect noticeably the classification into regu- 
lar/irregular types either. 



4.3 Classification according to star formation 
activity 

A typical problem of studies based on red galaxy samples is 
to disentangl e dust-reddened star-forming galaxies from qui- 
escent ones dPozzetti fc Mannuccill2000l ; iDaddi et all |2004 
iFaber et al1l2007l ; iPatton et al.ll201ll ). because both galaxy 
populations are indistinguishable us ing only broad-b and 
photometry at wavelengths < 10 fim i|Stern et al.ll2006l ). In 
particular, the rest-frame U — B color cannot differentiate 
between both red galaxy types adequately at z > 0.8 (F07). 
Therefore, different selection techniques based on color in- 
dices, mid-IR data, or SED fitting have been de veloped to 



isolat e both population s in red galaxy samples (iLin et al 



19991 ; IWolf et al. | 120031 ; IZucca et all l200rj ; IWilliams et a 



2009). 



The mid-IR emission and the S FR of a galax ; 
are known to be tightly correlated (|Kennicuttl 1 199 
IPerez-Gonzalez et all 120051 ). The higher sensitivities and 
spatial resolutions achieved by IR instruments in the 
last years have allowed the development of SFR indica- 
tors base d on the emission o f a ga l axy in a single mid- 
IR b and (ICalzetti et all 120051 . 120071 ; IPerez-Gonzalez et all 
2006). In particular, the 24 jim band of the Multiband Imag- 
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Figure 9. SFR cuts selected to define the red galaxies with en- 
hanced SFRs compared to the average SFR of the galaxy pop- 
ulation, as a function of redshift. Red squares: Redshift evolu- 
tio n of the SFR of ga laxies with log(.M*/M0) > 11.1 derived 
bv lKarim et alj {2011). Yellow circles: Redshift evolution of the 
SFR of galaxies with 10.8 < log(_M*/M Q ) < 11.1 derived by the 
same authors. Grey stars: Galaxy data from Perez- Gonzalez et al. 
(2005) with F(24fj,m) > 85 /xJ. Black stars: Galaxy data from the 
same authors with F(24fim) < 85 fj,]. This limiting flux in 24 (im 
naturally isolates galaxies with high (enhanced) SFR compared to 
the average value of the whole galaxy population at each redshift 
from those with low SFR, for the mass range considered here. 



ing Photometer in the Spitzer Space Telescope (MIPS) is 
found to be a good tracer of the infrared emission com- 
ing from the dust heated by star-forming stellar populations 
i|Alonso-Herrero et alll2006h . 

In this study, we have identified galaxies with notice- 
able star formation compared to the average SFR exhib- 
ited by the galaxies with similar masses at the same red- 
shift, because this is an evidence of that mechanisms dif- 
ferent to passive evolution are triggering it (such as tidal 
interactions, mergers, gas infall, or stripping). The SFR of 
a galaxy changes noticeably with redshift due to the natu- 
ral evolution of their stellar populations. The specific SFR 
has decayed with cosmic time as ~ (1 + z) n since z — 3, 
being n = 4.3 for all galax ies and n = 3.5 for star-forming 
sources (|Karim et al.ll201l[ ). This means that we must have 
into account the intrinsic rise of SFR with redshift to de- 
fine when a galaxy is forming stars more efficiently than 
the average of the whole galaxy population at each redshift. 
Therefore, we have used Spitzer/MIPS 24 fim data to dis- 
criminate between red galaxies with enhanced SFRs from 
those that have lower SFRs compared to the average SFR 
of the whole galaxy population at each redshift (see also 
I Williams et~alfl2009l '). A limiting flux in 24 /im of ~ 60 /xJy 
corres ponds to the 5-cr d etection level on our MIPS 24 /jm 
data (|Barro et al.ll2011ah . 

In Fig. 14. 2. 21 we plot the 2-evolution of the average SFR 
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Table 2. Percentages of AGN contamination in the subsample of 
HSF red galaxies by morphological types. 
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a Columns labelled as "Pure" represent the contribution of pure 
AGNs (Type 1 and 2) to each HSF galaxy type (i.e., without star- 
burst connection). 

b Columns labelled as "Total" provide the total contribution of AGNs 
to each type, including pure AGNs (Type 1 and 2) and mixed 
types (starburst-dominated AGNs, starburst-contaminatcd AGNs, 
and normal galaxies hosting AGNs). 

of galaxy populations with different mass ranges overlap- 
pi ng with ours ( which is \og(M* /M©) > 10.7), as derived 
bv lKarim et al.l l|201ll ) from a deep 3.6 /im-selected sample 
in the COSMOS field. We have overplot t ed th e location of 
the galaxies from I Perez- Gonzalez et ail (|2005T l in the dia- 
gram, differentiating those with emission fluxes in MIPS 
24 fim above 85 /iJ from those below it. Note that this lim- 
iting flux follows pretty well the redshift evolution of the 
average SFR of the whole galaxy population for our mass 
range, being a straightforward way for distinguishing galax- 
ies with enhanced SFRs from those with low SFRs, com- 
pared to the average value of the galaxy population at each 
redshift. Therefore, we have considered red galaxies with a 
24 /im MIPS flux above 85 ^iJy as systems with enhanced 
star formation compared to the global galaxy population 
at the same redshift (High Star-Forming galaxies or HSF, 
hereafter), whereas the galaxies with a 24 fim flux below 
this limit are not substantially star-forming compared to it 
(so, they are Low Star- Forming ones or LSF, henceforth). 

Active Galactic Nuclei (AGN) emission can coexist with 
star formation in galaxies, i.e., both phenomena are not mu- 
tually exclusive. Moreover, AGN activity is expected to be 
induce d and enhanced in gas-rich mergers (|Hopk ins et alj 
200830), so we should expect to find AGNs in our red galaxy 
sample, especially associated to our irregular class. Galaxies 
with an AGN are expected to be bright in the mid-to-far 
IR (and, in particular, in 24 /im) due to the emission of 
th e hot dust surrounding the central engine (see references 
in lPerez-Gonzalez et al.| [2005). Therefore, pure or dominant 
AGN emitters can be confused with HSF galaxies according 
to the criterion used in this study to identify HSF galaxies. 

Therefore, we have cross-correlated our red galaxy sam- 
ple !whh_th^_AG^£atalogue_over the Extended Groth Strip 
bv lRamos Almeida et al.Tl|2009l ) to estimate the contamina- 
tion of our HSF class by pure or dominant AGNs. These 
authors use a reliable technique of classification of active 
galactic nuclei (AGNs) based on fits of the UV-to-FIR SEDs 



of the galaxies with a complete set of AGN and starburst 
galaxy templates, in such a way that they are capable of 
distinguishing between pure and host-dominated AGNs up 
to z ~ 2 . The contribution of AGNs of any type and of 
pure AGNs (i.e., with no contribution of starbursts) to our 
HSF red galaxy sample at each redshift bin is tabulated in 
Tabled We also list the contribution of AGNs to HSF red 
galaxies at each redshift bin by morphological types. 

The results in Table[2]show that the contamination by 
any type of AGNs to our HSFs is negligible (below 13% at all 
redshifts), the contribution of pure AGNs being even smaller 
(< 7% at all redshifts). Therefore, we can conclude that 
AGN contamination does not affect our results concerning 
the classification into HSF galaxies significantly. 



5 ERRORS AND UNCERTAINTIES 

The errors in the estimates of the number density of each 
red galaxy type are derived considering the following sources 
of errors: statistical and classification errors, photometric 
redshift errors, and cosmic variance uncertainties. 

Statistical counting errors have been estimated for the 
number densities derived b y each classifie r for each galaxy 
type and at each redshift (|Gehreldll986h . The final num- 
ber density of each morphological type at each redshift bin 
is estimated as the mean of the number densities resulting 
from the three independent classifiers. We have estimated 
the error of this mean as the quadratic propagation of th e 
statistical error of each classifier (consult I Jogee et al.ll2009h . 

Concerning to redshift errors, we have used the red- 
shift estimates from the Rainbow Extragalactic Database 
|Barro et al.ll2011bl) . As commented in uncertainties in 
the redshifts are Az/(1 + z) < 0.03 for the whole redshift 
interval under consideration, for both bright and faint red 
sources. 

We have included estimates of the uncertainties in- 
troduced in the number densities of each red galaxy type 
by the redshift errors. These estimates have been derived 
using Monte-Carlo simulations. We have made 100 simu- 
lated catalogues of the (mass-limited) red galaxy sample, 
adopting a photometric redshift value for each source be- 
tween [z phot - A(z phot ), z phot + A(z phot )] at random, being 
Zphot the nominal photometric redshift of the source and 
A(zphot) the typical dispersion of the photometric redshifts 
compared to the spectroscopic ones. This dispersion is esti- 
mated as: A(zphot) =< A(z)/(1 + z) > -(1 + Zphot), where 
< A(z)/(1 + z) > is the average value obtained for this 
normalized dispersion at the redshift bin of the galaxy (see 
Fig-Hi. Then, we have obtained the number densities corre- 
sponding to each simulated catalogue for each galaxy type 
and at each redshift bin, accounting for the different red- 
shifts of each catalogue. The dispersion of the 100 values 
obtained for the number density at each redshift bin and 
galaxy type represents an estimate of the error associated 
to the photometric redshift uncertainties. 

Statistics of massive red galaxies can be dramatically 
affected by cosmic variance due to their high clustering 
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|Somerville et al.ll2004). We have estimated cosmic variance 
using the model bv lMoster et al.l |201ll 'l. which provides es- 
timates of cosmic variance for a given galaxy population us- 
ing predictions from cold dark matter theory and the galaxy 
bias. They have developed a simple recipe to compute cosmic 
variance for a survey as a function of the angular dimensions 
of the field and its geometry, the mean redshift and the width 
of the considered redshift interval, and the stellar mass of the 
galaxy population. We have considered the geometry and an- 
gular dimensions of our field, as well as the different redshift 
bins analysed in each case to estimate the cosmic variance. 
iMoster et al.l software provides these estimates in two mass 
ranges overlapping with ours: 10.5 < log(.M*/M0) < 11 
and 11 < log(A4, /Mq) < 11.5. Therefore, we have consid- 
ered the mean cosmic variance of both mass ranges as a rep- 
resentative value of the cosmic variance of our mass-limited 
sample at each redshift bin. Cosmic variance depends on the 
redshift. At 0.3 < z < 1.5 and for our mass range, cosmic 
variance decreases from 0.3 < z < 0.5 to 0.7 < z < 0.9 
approximately, and then starts rising again towards higher 
redshifts. This behaviour is observed for all used redshift 
bins (Az = 0.1, 0.2, and 0.4). For Az = 0.1, the root cosmic 
variance acquires minimum and maximum values equal to 
34% and 40%, respectively. For Az = 0.2, it changes between 
24% and 28%, and between 17% and 19% for Az = 0.4. All 
our results include these uncertainties quadratically added 
to the other error sources (statistical and classification er- 
rors, and redshift uncertainties). 



The predi ctions of the cosmic variance obtained with 
the model by IMoster et al.l ()201 lh are in good agreement 
with the rough estimates that can be derived from panel 
a of Fig. [5] At 0.7 < z < 1.1, the dispersion of the dif- 
ferent studies is ~ 30%, most of which can be attributed 
to cosmic variance. This value is quite similar to the esti- 
mate obtained with IMoster et al.l model at these redshifts 
for Az = 0.2 (~ 27%). At lower and higher redshifts, the 
dispersion among different authors may arise also in com- 
pleteness problems and higher observational errors, so a di- 
rect comparison cannot be done. 



The redshift errors have been added quadratically to the 
statistical and classification errors and to the uncertainties 
associated to cosmic variance to obtain the final error of 
the number density at each z interval and for each galaxy 
type. In general, for regular red galaxies, the statistical and 
classification errors contribute to ~ 30% to the total error 
of their number density for both low and high redshifts, the 
redshift errors represent ~ 10% of the total error, and the 
cosmic variance contributes to ~ 60% of the total error at 
all redshifts. For irregular galaxies, the contribution of the 
statistical and morphological classification errors represents 
~ 35% at all redshifts, the redshift errors are ~ 25% (being 
~ 30% at low z, and decreasing down to ~ 23% at high z), 
and the cosmic variance contributes to ~ 37% (being ~ 35% 
at low z, and rising to ~ 40% at high z). 



6 OBSERVATIONAL TESTS TO THE 

HIERARCHICAL ORIGIN OF MASSIVE 
E-SO'S 

As commented in JTJ no observational evidence on the exis- 
tence of a evolutionary link between major mergers and the 
rise of the present-day massive E-SO's exists up to the date. 
In order to test the major-merger origin of massive E-SO's 
observationally, we must be capable of providing evidence 
supporting or rejecting the existence of such a link. How- 
ever, the transitory stages and end products of these evo- 
lutionary tracks coexist at each redshift, making this task 
difficult. Here we define three tests to observational data 
based on the predictions of hierarchical models to check if 
data are coherent with the existence of this link, that can 
be done thanks to the classification performed in this study. 

6.1 Equivalence between our galaxy types and the 
different evolutionary stages of red galaxies 

The morphological and structural properties of a galaxy can 
be combined with information about its star formation en- 
hancement to establish correspondences between its galaxy 
type in our classification and the evolutionary stage of the 
galaxy. Once a galaxy is identified as a major merger on 
the basis of its noticeable structural distortion, both its star 
formation level and its global morphology may be pointing 
either to the gas-content of its progenitors or to the merger 
phase (|Cox et al.ll2008l ; lLotz et al.ll2008bl . l2010al fbh . 

Galaxies involved in a gas-rich major merger are ex- 
pected to be disk-dominated during the merging-nuclei 
phase. During the initial moments of this phase, these 
mergers experience strong starbursts that last less than 
~ 0.5 Gyr. Therefore, depending on the evolutionary stage 
of these mergers, they might appear either as HSF or 
LSF. Post-merger stages of these events will still exhibit 
noticeable distortion, but its morphology is expected to 
be spheroidal-dominated. Depending on the efficiency of 
the star formation quenching in the merger, late phases 
of gas-rich major mergers may be HSF or LSF. Addi- 
tionally, gas-poor mergers are expected to be irregular, 
spheroidal-dominated, and mostly LSF since their merging- 
nuclei phase. Normal disks and E-SO's exhibit a regular disk 
or spheroidal-dominated morphology, respectively. Most E- 
SO's are quiescent, meaning that they must be LSF. 

Our classification thus allows us to distinguish among 
these evolutionary stages of red galaxies just attending to 
their types, except between gas-poor mergers and post- 
merger stages of wet mergers. In Table[3] we list the cor- 
respondence between a galaxy evolutionary stage and its 
characteristics in our classification. 

6.2 Hierarchical evolutionary paths among red 
galaxy types 

Here, we describe three tests based on the expectations of 
hierarchical models of galaxy formation to check whether 
data is coherent with the existence of an evolutionary link 
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Table 3. Correspondence between different galaxy evolutionary stages and our red galaxy types. 



Galaxy Type 


Global morphology 


Structural distortion a 


Star formation 13 


(1) 


(2) 


(3) 


(4) 


E-SO 


Spheroid-dominated 


Regular 


LSF 


Spiral 


Disk-dominated 


Regular 


HSF/LSF 


Gas-rich major mergers (merging-nuclei phase) 


Disk-dominated 


Irregular 


HSF/LSF 


Gas-rich major mergers (post-merger phase) 


Spheroid-dominated 


Irregular 


HSF/LSF 


Gas-poor major mergers 


Spheroid-dominated 


Irregular 


LSF 



a Non-distorted morphologies are noted as regular, and merger-like distorted ones as irregulars. 

b HSF: objects with enhanced star formation compared to the average of the whole galaxy population; LSF: objects with star formation 
lower than this average SFR. 

c Only major mergers already merged into one body are traced by this study. Galaxy pairs are considered as two still-independent galaxies. 



between major mergers and the E-SO's appearing on the Red 
Sequence since z ~ 1.5 or not. 

F07 proposed a hierarchical mixed evolutionary sce- 
nario for explaining the observed mass migration from the 
massive end of the Blue Galaxy Cloud to that of the Red 
Sequence since z ~ 2, in which "quenched galaxies enter the 
Red Sequence via gas-rich mergers" , and can be "followed by 
a limited number of gas-poor, stellar m ergers along the se- 
quence " . The semi-analytical model by lEliche-Moral et"al] 
l|2010al . EM10 hereafter) proved the feasibility of the F07 
scenario for explaining the buildup of E-SO's that end up 
with M* > 10 Mq at z = 0, just accounting for the effects 
of the m ajor mergers strictly report ed by observations since 
z ~ 1.2 (|L6pez-Saniuan et al.r 2009). This model reproduces 
the observed evolution of the massive end of the galaxy lumi- 
nosity function by color and morphological types. The evo- 
lutionary track described by F07 appears naturally in the 
model, as it considers the relative contribution o f gas-poor 
and g as-rich mergers at each redshift reported bv lLin et alj 
(|2008h and their different effects on galaxy evolution. 

The advantage of this model is that its predictions are 
in excellent agreement with cosmological hierarchical mod- 
els (despite being based on observational major merger frac- 
tions) , reproducing observational data at the same time (see 
EM10; lEliche- Moral et al.l l2010bl ). Based on these predic- 
tions, we have defined some tests that observational data 
must fulfill if most massive E-SO's have really derived from 
major mergers occurred at relatively late epochs in the cos- 
mic history. These predictions are the following ones: 

(i) Most present-day E-SO's with M» > lO n M are the 
result of at least one gas-rich major merger that place them 
on the Red Sequence since z ~ 1.2. 

(ii) In addition, ~ 75% of the remnants resulting from 
these gas-rich events have been involved in a subsequent 
gas-poor major merger, occurred quite immediately. The re- 
maining ~ 25% have thus continued their evolution towards 
an E-SO passing through a quiet post-merger phase. 

(iii) The bulk of these major mergers are at intermediate- 
to-late stages during the ~ 2 Gyr period elapsed at 0.7 < z < 
1.2, which means that the gas-rich ones must have started 
at 1 < z < 1.5, accouting for the typical timescales of these 



events. The gas-poor ones are later, but must take place 
earlier than z ~ 0.7 in any case, as the resulting massive 
E-SO's have been in place since that epoch. 

Therefore, according to the EM10 model, the appear- 
ance of the bulk of massive E-SO's takes place at 0.7 < z < 
1.2, and nearly all have evolved according to the following 
path: 



Gas-rich major merger 



Post-merger stage 

[ ~ 25% 1 



Gas-poor major merger 

[ ~ 75% 1 



Massive E-SOa 



(1) 



Note that this evolutionary track implies the existence of a 
nearly 1:1:1 numerical relation at 0.7 < z < 1.2 between gas- 
rich major mergers at merging-nuclei stages, gas-poor events 
and post-mergers stages of gas-rich ones, and the massive 
E-SO's assembled at those epochs. Accounting for the corre- 
spondence of these evolutionary stages and the red galaxy 
types denned in this study (see Table[3J), the evolutionary 
path schematized in eq.[T]can be re- written as follows: 



ID 

(HSF/LSF) 



IS (HSF/LSF) 

[ ~ 25% 1 



IS (mostly LSF) 

[ ~ 75% 1 



4 



RS (mostly LSF) 



(2) 



Then, the previous 1:1:1 relation between the different stages 
of these galaxies in their evolution towards the Red Sequence 
can be translated into an equivalent 1:1:1 relation between 
the following red galaxy types, according to Table[3] 



ID -> IS 



RS. 



(3) 
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If most massive E-SO's have really evolved according to the 
hierarchical scenario proposed by the EM10 model, then the 
data must fulfill this evolutionary path among massive red 
galaxy types at 0.7 < z < 1.2. This imposes the following 
three observational tests or constraints to observations: 

1. The accumulated number density of ID's (gas-rich ma- 
jor mergers) on our red sample since z ~ 1.2 down to a lower 
redshift z must reproduce the net numerical increment of 
RS's (E-SO's) observed between the two redshifts. 

2. An analogous relation must be fulfilled by the accumu- 
lated number density of IS's (gas-poor major mergers and 
post-merger stages of gas-rich events). 

3. The bulk of RS's (E-SO's) with stellar masses M. > 
1O 11 M0 at 2 = must have been definitely built up during 
the ~ 2.2 Gyr time period elapsed at 0.7 < z < 1.2. 

Our data is complete for galaxy masses M» > 5 x 
1O 1O M at each redshift. This means that we can ensure 
that we are in a position to trace back in time the potential 
progenitors of the present-day E-SO's with M» > 10 11 Mq at 
z ~ that could have merged to create them during the last 
~ 9 Gyr (see 

Therefore, we have estimated the cumulative distribu- 
tion of ID's and of IS's, and we have compared them with 
the redshift evolution of the number density of RS's since 
z ~ 1.5. In the case that major mergers have not driven the 
assembly of the massive E-SO's as proposed by the EM10 
model, the previous three tests must fail. On the contrary, 
if these three distributions agree pretty well, these tests will 
support strongly the existence of an evolutionary link be- 
tween major mergers and the appearance of massive E-SO's, 
as expected by hierarchical scenarios of galaxy formation. 
The results of these tests are presented in i]7.3l 

We must remark that the EM10 model exclusively 
quantifies the effects of major mergers on galaxy evolu- 
tion at z < 1.2. Hence, it does not discard the contri- 
bution of different evolutionary processes to the definitive 
assembly of massive E-SO's, although it predicts that it 
must have been low, as most of their number density evo- 
lution can be explained just accounting for the effects of 
the major mergers. This seems to be confirmed by obser- 
vations, as other evolutionary mechanisms (such as minor 
mergers, ram pressure stripping, or bars) seem to have 
been significant for the formation of the Red Sequence 
only for low and intermediate masses , but not fo r high 



2007 



masses (Kormcndy & Kcnnicutt 
Domuigucz-Pa lmero fc Balcellsl 120 08; Bu ndy 
Simard et ail 120091 ; ICameron et all |2010| ; iKavirai et al.l 



2009 



201 if ). Moreover, the model does not exclude disk rebuilding 



after the major merger either. On the contrary, it is proba- 
bly required for giving rise t o a SO instead of an elliptic al, 
as indicated by o bservations ( Ha mmer et al.ll2005l . l2009al lbt 
I Yang et al.ll2009t ). 

Moreover, the EM10 model assumes that intermediate- 
to-late stages of major mergers are red and will produce 
an E-S0, on the basis of many observation al and compu- 
tational studies (see references in EM 10 and Cimatti et alj 



These assumptions are crucial for the model, as they are nec- 
essary to reproduce the redshift evolution of the luminosity 
functions selected by color and morphological type. There- 
fore, testing if our data is coherent with the existence of 
an evolutionary link between the advanced stages of major 
mergers in our red sample and the definitive buildup of mas- 
sive E-SO's, we are also indirectly testing these assumptions 
of the EM10 model. 



6.3 Observational considerations for the tests 

According to iBell et all (|2006l ). the average number density 
of major mergers at a given redshift centered at z, n m (z), 
is related to the number density of the major mergers de- 
tected at certain intermediate phase of the encounter in that 
redshift bin, < n ul>p i lasc (z) >, as follows: 



n m (z) =< n m ,p haso (,z) > 



t([zi,z 2 ]) 



Tdct 



(4) 



l2002l : ISchawinski et al 1l2007l , l2010l : IChilingarian et al.ll2010h 



with t( [21,22]) being the time elapsed in the redshift bin, 
and Tdot representing the detectability time of the interme- 
diate merger stage under consideration. We have used this 
equation in our tests. 

We find that the number densities of ID's and IS's re- 
main quite constant with redshift (see iJTJ). As ID's and IS's 
correspond to intermediate merger stages, this means that 
the major merger rate must evolve smoothly with redshift, in 
good agreement with observational estim a tes of merger rate s 
l|Lin et all 120081 ; iBrammer et all l2011bl ; lLotz et alj l201ll ). 
This also indicates that the net flux of irregular galaxies 
appearing on the Red Sequence or at nearby locations on 
the Green Valley (i.e., the number of red irregulars created 
and destroyed per unit time in the red sample) must have 
been nearly constant at 0.3 < 2 < 1.5 for our mass range. 
Together with the evolutionary track stated in eq. [3] this 
implies that the accumulated number of ID's observed in a 
redshift bin must be equal to the accumulated number of 
IS's in the same bin, and also to the number of RS's as- 
sembled at the end of the redshift bin. As we also find that 
the number of RS's is negligible at 2 = 1.5 (see 33, it is 
also justified to start accumulating ID's and IS's since this 
redshift down to 2 ~ 0.3. 

In order to estim ate n m (g), we have adopted the merger 
timescales derived bv lLotz et all (|2010al ). These authors re- 
port typical detectability timescales through CAS indices 
of simulated major mergers, as a function of the baryonic 
gas fraction and the mass ratio of the encounter. For gas 
fractions represent ative of gas-rich merg ers up to 2 ~ 1 
(~ 30 - 50%, see H ammer et all l2009af ) and mass ratios 
typical of major mergers (1:1 to 1:4), they find an aver- 
age detectability timescale Tdet, gas-rich ~ 1.0 ± 0.3 Gyr of 
merging-nuclei phases. For gas fractions representative of 
gas-poor major mergers (< 30%), we estimated from these 
author and from Ivan Dokkuml l|2005h an average timescale 
of Tdet, gas-poor ~ 0.5±0.2 Gyr of their merging-nuclei phase. 
The timescales for post-merger stages are similar to those of 
their me rging-nuclei phases for both gas-rich and ga s-poor 
mergers l|Schweizer fc Seitzerll2007l ; lLotz et al.ll2010al ). The 
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errors in these timescales account for small chang es of these 
repres entative values with redshift, accordingly to lLin et al.l 

|2oid y 



7 RESULTS 

In TableUwe show the number densities of red galaxies with 
M* > 5 x 1O 1O M per redshift bins at 0.3 < z < 1.5 obtained 
from our data, according to the classification resulting from 
the morphological/distortion level of the galaxy f N4.HI and 
from its star formation level ( i]4.3|) . Therefore, the galax- 
ies are classified as regulars or irregulars (according to their 
structural distortion level), as spheroid/disk-dominated (ac- 
cording to their global morphology), and as HSF/LSF (ac- 
cording to their star formation enhancement). Errors in the 
results listed in this table and in all the figures of the present 
section account for the statistical and classification errors, 
the redshift errors, and the cosmic variance uncertainties, as 
described in iJS] 



7.1 Morphological evolution of red galaxies since 

z ~ 1.5 

Figure[TU] shows the redshift evolution of the fraction and 
the comoving number density of massive red galaxies by 
galaxy types considering their global morphology and struc- 
tural distortion level. The fractions and number densities of 
red regular galaxies rise with cosmic time, with the fraction 
of spheroids being higher than that of disks at the same red- 
shift. Therefore, we are tracing the progressive settlement of 
regular galaxies on the massive end of the Red Sequence at 
0.3 < z < 1.5, made primarily of spheroids but also contain- 
ing some regular disks (mainly at z < 0.7). We also find sig- 
nificant populations of irregular galaxies at all redshifts up 
to z ~ 1.5, with number densities > 10 -4 Mpc -3 . Although 
the fractions of red irregular spheroids and disks decrease 
with cosmic time, their densities remain quite constant at 
0.3 < z < 1.5. This fact clearly points to the transitory na- 
ture of red irregular galaxies at any redshift, as indicated by 
previous studies l|Lotz et al.ll2008al ). 

Figure[TT] shows the redshift evolution of the number 
density and fractions of red regular galaxies (RS + RD) and 
irregulars (IS + ID) at 0.3 < z < 1.5, for the three wide 
redshift bins under consideration. The fraction of regular 
galaxies increases by a factor of ~ 6 since z ~ 1.3 down 
to z ~ 0.5, whereas the fraction of irregulars decreases by 
the same amount (top panel in the figure). Moreover, the 
number density of regulars has risen by a factor of ~ 12 
during this time period (bottom panel in the figure), while 
that of the irregulars keeps constant until z ~ 0.9, decreasing 
by a factor of ~ 3 at z < 0.6 ( in excellent agreement with 
the results bv lllbert et~al]|201oh . 

At z ~ 1.3, irregular galaxies (i.e., major mergers) rep- 
resent nearly 80% of the red galaxy population. Their frac- 
tion decreases down to ~ 50% at z ~ 0.9, and reaches 
< 15% at z ~ 0.5. Considering that ~ 70% of these ir- 
regular galaxies are dust-reddened due to intense star for- 
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Figure 11. Redshift evolution of the number density and fraction 
of all regular and irregular (i.e., major mergers) red galaxies with 
M, > 5 X 10 10 M Q at 0.3 < z < 1.5, in the three wide redshift 
bins under consideration (0.3 < z < 0.7, 0.7 < z < 1.1, and 
1.1 < z < 1.5). Top panel: Redshift evolution of the fraction 
of red massive regular and irregular galaxies with respect to the 
total red galaxy population. Bottom panel: Redshift evolution of 
the number density of regular and irregular massive red galaxies 
in the same redshift bins. 



mation (see H7.2p . these results are in good agreement with 
previous estimates of the fraction of dust-reddened star- 
forming galaxies on the Red Sequence at dif f erent redshifts 
dCimatti et al.ll2002l; lYan fc Thompsonll2003!; jGilbank et al ' 
| 2003 |; i Bell et alj|2004al; iMoustakas et alj|2004 IWeiner et al 
20051 ; Icirasuolo et alj|2007l ; iFranzetti et al.ll2007l ; IZhu et al 



We have derived the redshift evolution of the number 
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Table 4. Comoving number densities of the different red galaxy types in units of xlO - 4 Mpc~ 3 . 
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a Net number of galaxies in the final sample in a sky area of ~ 155arcmin 2 . 

b Galaxy types according to morphology: Regular Spheroid-dominated (RS), Regular Disk-dominated (RD), Irregular Spheroid-dominated 
(IS), Irregular Disk-dominated (ID). 

c Galaxy types according to star formation activity: High Star-Forming (HSF), Low Star-Forming (LSF). 
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Figure 10. Redshift evolution of the number density and fraction of red galaxies with M* > 5 X 10 Mq, according to their morphology 
and structural distortion level, for the three wide redshift bins considered in the study (0.3 < z < 0.7, 0.7 < z < 1.1, and 1.1 < z < 1.5). 
Top panels: Redshift evolution of the fraction of each type with respect to the whole red galaxy population at each redshift bin. Bottom 
panels: Redshift evolution of the number density of each morphological type in the same redshift bins. 



density of regular and irregular massive red galaxies in nar- 
rower redshift bins, to delimit more accurately the epoch at 
which the bulk of the red regular galaxies appears into the 
cosmic scenario (see Fig. ll2p . The large uncertainties of our 
results in these narrow redshift bins prevent us of deriving 



quantitative conclusions based on this figure. However, we 
can conclude from it that the number density of red irreg- 
ulars (i.e., major mergers) has remained nearly constant at 
0." 7 z J; 1-5 within errors. This indicates that these sys- 
tems must have been forming and disappearing at similar 
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Figure 12. Redshift evolution of the number density of red reg- 
ulars and irregulars with M, > 5 X 10 10 Mq at 0.3 < z < 1.5, 
using narrow redshift bins of Az = 0.1. 



rates during this time period. Therefore, the decrease of the 
fraction of irregular systems with cosmic time is exclusively 
relative, i.e., it is just due to the rise of the population of reg- 
ular galaxies, but not to a net drop in their number density 
down to z ~ 0.6. The figure also shows that the redshift- 
evolution of the number densities of regulars and irregulars 
on the Red Sequence cross at z = 0.9 for M* > 5 x 10 10 M Q . 

As commented in JJT] the rise of the number density 
of red regular types with cosmic time and the constancy of 
that of irregular ones has been interpreted as a sign pointing 
to the conversion of irregulars into regulars with time. We 
provide observational evidence supporting the existence of 
this evolutionary link in Q7.3\ 



7.2 Star formation activity according to red 
galaxy types since z ~ 1.5 

In Fig.[T2]we show the redshift evolution of the number den- 
sity of massive red galaxies for each morphological type de- 
fined in this study, attending to its star formation activ- 
ity. We remark that HSF's are galaxies that show enhanced 
SFRs compared to the average SFR of the galaxy popula- 
tion at each redshift f i|4.3[) . All types, except RS's, host a 
significant number of HSF galaxies, with percentages vary- 
ing depending on type and redshift. Red regular spheroids 
are the galaxy types hosting the lowest fractions of HSF sys- 
tems at all redshifts since z ~ 1.5, as expected. Curiously, 
red RD's exhibit a noticeable increase of the HSF systems 
fraction at z < 0.7 (panel b). Irregular types harbour en- 
hanced SFRs typically (both spheroidal and disk systems), 



coherently with their merger-related nature (panels c and d 
in the figure) . The percentage of HSF objects in these types 
has not changed at 0.3 < z < 1.5 within errors, the fraction 
of HSF objects in ID's being much higher than in IS's at 
all redshifts (~ 80% for ID's and ~ 50% for IS's). This is 
normal if we consider that the former ones correspond to 
intermediate stages of gas-rich mergers, while the later ones 
group gas-poor major mergers and post-merger phases of 
the gas-rich ones. 

Figure[T3l implies that, at 0.7 < z < 1.5, the enhanced 
star formation in red galaxies is mostly hosted by major 
mergers; mostly by gas-rich merger remnants at intermedi- 
ate phases (ID's, which represents ~ 30% of all red galaxies 
at these redshifts, see Fig. [TO} and by ~ 50% of spheroidal 
remnants (IS's, which represent ~ 30%). At z < 0.7, en- 
hanced SFRs can be found still in most red major mergers 
(which represent ~ 25% of the whole red population, see 
Fig.H0) and in nearly all red RD's (~ 15% of it). RS's are 
the types that host the lowest levels of SFR enhancement 
at all redshifts, hosting ~ 20% at most at z < 0.7 (despite 
being ~ 60% of the whole red population at z < 0.7). 

Considering that our detectability timescale for gas-rich 
mergers is r ~ 1 Gyr ( >i6.3|l , the merger rate derived from our 
data for ID's at 0.3 < z < 0.7 is ~ 2.5 x 10~ 4 Gyr" 1 Mpc" 3 
(considering their number density from Fig.ll3[). which is an 
est imate that i s in g ood agreement with the one reported 
bv lChou et ail (|201l( ) for gas- rich mergers using a different 
methodology (see their Figure 6) . This strongly supports the 
identification of our ID's with intermediate stages of gas-rich 
major mergers. 



7.3 Tests to the definitive buildup of massive 
E-SO's through major mergers 

This section presents the results of the three observational 
tests commented in H7.3I to test the hierarchical formation 
scenario of massive E-SO's. 



7.3.1 Test 1: Cumulative redshift distributions of red 
irregular disks since z ~ 1.5 

The first panel of Fig. [14] shows the result of the first test 
proposed in H6.2I We plot the cumulative number density 
distributions of red ID's (intermediate stages of gas-rich ma- 
jor mergers, see Table[3]l since z ~ 1.5 down to redshift z as 
obtained from eq.[3] and compare it with the redshift evo- 
lution of the number density of red RS's (E-SO's, see the 
same table). The detectability timescale considered for the 
ID's is that commented in H6.3I The filled area around the 
cumulative data points accounts for the uncertainties in this 
timescale. The number densities at the lowest redshift bin 
are probably affected by volume and selection effects as sug- 
gested by this figure and as indicated in £13.31 so we will 
override this redshift bin in all our results and comments 
henceforth. 

The cumulative distribution of ID's reproduces pretty 
well both the redshift evolution of all RS's (E-SO's) and that 
of exclusively the LSF RS's (quiescent E-SO's) at 0.6 < z < 
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Figure 13. Rcdshift evolution of the comoving number density of red galaxies with M* > 5 X 10 10 Mq, according to their morphological, 
structural, and star formation activity properties. Solid lines: Number densities for LSF galaxies of each type. Dashed lines: Number 
densities of HSF galaxies of each morphological type. Panel a): Regular spheroids (RS's). Panel b): Regular disks (RD's). Panel c): 
Irregular spheroids (IS's). Panel d): Irregular disks (ID's). 
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Figure 14. Observational tests to the data of red galaxies for testing the hierarchical buildup of massive E-SO's, on the basis of the 
expectations of the EM10 model. Panel a) Test 1: Cumulative redshift distribution of red ID's (intcrmcdiate-to-late stages of gas-rich 
major mergers) since z = 1.5 down to z = 0.3, compared to the buildup of massive RS's (E-SO's) during the same time period. Panel 
b) Test 2: Cumulative redshift distribution of red IS's (intermediate-to-late stages of gas-poor major mergers and late stages of gas-rich 
ones) since z = 1.5 down to z = 0.3, compared to the buildup of massive RS's (E-SO's) during the same time period. Panel c) Test 
3: Comparison of the redshift evolution of the number density of RS's (E-SO's) for two different mass ranges (M» > 5 x 1O 1O M and 
M * ^ 10 1:l Mq) and the predictions of the EM10 model for E-SO's that end up with M» > 10 11 M Q at z = 0. Vertical solid line: Epoch 
at which the E-SO's with M, > 10 n M Q at z = can be considered as definitively asse mbled according to the EM10 model. Blue circle: 
Local number density of red galaxies with M» > IO^Mq derived by iMadgwick et al] ll2003h from 2dF survey. The filled areas around 
the cumulative data points in panels a and b correspond to the uncertainties in the merger detectability timescales, while in panel c 
indicates the model uncertainties. The lowest redshift bin probably suffers from incompleteness. 



1.5 within errors (most E-SO's at these redshifts are LSF's, 
see Fig, ll3[) . This is compatible with that the bulk of red 
gas-rich major mergers at intermediate stages of their evo- 
lution at these redshifts evolve into the E-SO's that start to 
populate the Red Sequence at later cosmic times. Therefore, 
panel a in Fig.[l4]supports the existence of the following evo- 
lutionary link among these two massive red galaxy types at 
0.6 < z < 1.5: 

ID -> RS. (5) 
Panel a of Fig.fJJ] provides observational support to 



many aspects concerning the buildup of massive red galaxies 
expected by hierarchical models of galaxy formation: 

(i) At 0.6 < z < 1.5 all irregular disks populating the 
massive end of the Red Sequence and nearby locations have 
evolved into the regular spheroids that appear on it some 
time later. 

(ii) This implies that the bulk of red massive E-SO's at 
z ~ 0.6 derive from disk galaxies that have migrated from 
the Blue Cloud to the Red Sequence during the time period 
enclosed at 0.6 < z < 1.2. 

(iii) It also supports that the dominant mechanism after 
this evolution has been major merging, as data is compatible 
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with the fact that most massive E-SO's that have appeared 
on the Red Sequence since z ~ 1.2 seem to have undergone 
a previous phase of gas-rich major merger (i.e., a previous 
evolutionary stage as a red ID). 

We remark that the major mergers detected at their 
phases as ID's by z ~ 1.2 in this study must have started ~ 
0.5-0.7 Gyr before, i.e., they were at their approaching (pre- 
merger) phases by z ~ 1.5. Therefore, the first generation 
of E-SO's that should result from these mergers according 
to the hierarchical scenario have appeared after z ~ 1.2. 
Accounting for this and for the fact that the EM10 model 
traces the evolution of massive E-SO's back-in-time up to 
z — 1.2, a comparison between our data and this model is 
feasible only up to this redshift. 

In conclusion, the panel a of Fig.[2] indicates that the 
data fulfill the first observational constraint imposed by the 
EM10 model ( i]6.2|l . supporting the evolutionary track indi- 
cated in eq. [T] (see also eqs(2] and [3} • 

7.3.2 Test 2: Cumulative redshift distributions of red 
irregular spheroids since z ~ 1.5 

Panel b of Fig.[l4] now compares the accumulated redshift 
distribution of all IS's (gas-poor mergers plus post-merger 
stages of gas-rich ones, see Table[3]) with the redshift evolu- 
tion of the number density of massive E-SO's (RS's) since 
z ~ 1.5. The detectability timescales considered for the 
IS's are those commented in £16.31 The filled area around 
the cumulative data points accounts for the uncertainties in 
this timescale. Again, the accumulated distribution of IS's 
reproduces quite nearly the buildup of massive E-SO's at 
0.6 < z < 1.5 within errors. Therefore, this figure sup- 
ports the existence of an additional evolutionary link at 
0.6 < z < 1.5 among the following red galaxy types: 

IS -> RS. (6) 

This panel of Fig. [14] also shows that data is compatible 
with the fact that most massive red E-SO's at z ~ 0.6 have 
also experienced a previous transitory phase as a red IS at 
0.6 < z < 1.5. The evolutionary links shown in eqs.[5] and [6] 
can be summarized into the following evolutionary track: 

ID -> IS -»■ RS, (7) 

which is the evolutionary path between red galaxy types 
expected by the hierarchical scenario of the EM10 model 
for the buildup of massive E-SO's (see eq.0). Therefore, the 
panel b of Fig. 1 141 shows that the data also fulfill the second 
constraint imposed by the hierarchical evolutionary scenario 
proposed in the EM10 model (~ t|6.2|) . supporting the evolu- 
tionary track indicated in eq. [T] 

7.3.3 Test 3: Epoch of definitive assembly of massive 
E-SO's 

In the last panel of Fig. 1141 we compare the redshift evolu- 
tion of the number density of E-SO's predicted by the EM10 
model with the one obtained from our data for two different 



mass selections: for our nominal one (M* > 5 x 1O 1O M0), 
and for an alternative one which is identical to our nominal 
one in all aspects, except in the mass range (M* > lO n M0). 

The EM10 model traces back in time the evolution of 
the E-SO's that have M* > IO^Mq at z = (MM- Ac- 
cording to the model, these galaxies have been mostly as- 
sembled through major mergers at 0.7 < z < 1.2, their 
buildup being frozen since then. Therefore, the model starts 
to trace the progenitors of these E-SO's at z > 0.7, which 
have masses lower by a factor of ~ 2 compared to the E-SO's 
resulting from the merger. Consequently, the model predic- 
tions on the number density of E-SO's at z > 0.9 can be 
compared with our results for M, > 5 x 1O 1O M0 (as both 
studies trace similar mass ranges globally). But at lower 
redshifts, the EM10 model traces E-SO's that already have 
M» > 1O 11 M0, so it is comparable to the results with a mass 
selection M* > 1O U M . 

As Fig. 1141 c shows, at z > 0.7, the model repro- 
duces much better the settlement of the RS's with M* > 
5 x 1O 1O M0, whereas at z < 0.7 it clearly follows the trend 
of RS's with M, > lO n M0 , as expected from the arguments 
given above. However, our data seem to have completeness 
problems at z < 0.5, so we cannot assure that the number 
density of E-SO's with M» > lO n M0 has remained constant 
since z ~ 0.7. 

Nevertheless, the number density of E-SO's with M« > 
10 n M© at z ~ 0.6 estimated with our data is quite sim- 
ilar to the one estim ated for these galaxies at z — 
l|Madgwick et alj|2003t ). as shown in panel c of Fig.lTl This 
means that the number density of these objects has re- 
mained nearly constant since z ~ 0.6, in good agreement 
with the predictions of the EM10 model. 

Although better data at z < 0.6 are required to directly 
confirm this result, our data are coherent with the fact that 
E-SO's with M, > IO^Mq have been definitively assembled 
since z ~ 0.6, supporting that the bulk of the assembly of 
red massive E-SO's has occurred during the ~ 2.2 Gyr period 
elapsed at 0.7 < z < 1.2, as predicted by the EM10 model. 
Therefore, we can conclude that our data fulfill the third 
constraint imposed by the hierarchical evolutionary scenario 
proposed in the EM10 model ( H6.2p , supporting again the 
evolutionary track indicated in eq. [1] 

Summarizing, the three observational tests to the hi- 
erarchical evolutionary framework proposed by the EM10 
model strongly supports this scenario, pointing to major 
mergers as the main mechanism for the buildup of massive 
E-SO's. 



8 DISCUSION 

In general, our results strongly support a late definitive for- 
mation redshift for massive E-S0 (z < 1.5), in agreement 
with hierarchi cal scenario s of ga laxy formation (see refer- 
ences m ©. Illbert et al.l (|2010T ) already claimed for the 
freezing in the assembly of massive E-SO's at z < 0.7, propos- 
ing that a sudden drop of the gas-rich major merger rate 
must take place at z < 0.8 to explain it. This fact was 
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proven to be observationally feasible by the EM10 model, 
just accountin g for the major merger fr actions reported by 
observations l|Eliche-Moral et alj|2010bh . The present study 
supports it observationally. 

Apart from the major mergers that seem to have been 
the main responsible mechanisms to place the massive E- 
SO's on the Red Sequence at 0.6 < z < 1.5, other evolu- 
tionary processes must have contributed to the evolution of 
massive E-SO's down to the present (and even coeval with 
them). But, as commented above, they must not have risen 
their masses noticeably or changed their morphology notice- 
ably, as E-SO's with log(M»/M0) > 11 seem to have been in 
place since z ~ 0.6 (see i|7.3[) . This is supported by the fact 
that many of these processes are observed to be relevant only 
for the evolution of galaxies with lower masses dBekki et al.l 
120021 ; iDriver et af]|l998l ; lAragon-Salamanca et alj|2006l in 
particular, the fading of stellar populations, see). 

We must also remark that our red galaxy sample does 
not trace the evolution of SO's observed in the clusters, 
which seems to have been relevant since z ~ 0.4-0.5 due to 
the environmental star-fo rmation quenching of the spirals 
that fall into the clusters jDesai et al.|[2 007; Poggiant i et al.l 
l2009bl ; ISimard et aT1l2009l ; IVulcani et al.ll2011i ). Cluster SO's 
have typical ma sses lower than those selected here (M« < 
5 x 1O 1O M , see iDressler et al. Hl999l ; lBedregal et alJ^OOrl . 
and hence, our results do not apply to them in general. More- 
over, some studies indicate that the fra ction of SOs in group s 
is similar to that of clusters at z < 0.5 l|Wilman et al.| [2009). 
Considering that most galaxies reside in groups (~ 70%, see 
iBerlind et aill2006l ; ICrook et al.ll2007l 1. this means that the 
majority of SOs in the Universe are located in groups (not 
in clusters). Nevertheless, note that this evolution in clus- 
ters does not contradict the EM10 model at all, because this 
model exclusively analyses the effects of the major mergers 
on galaxy evolution since z ~ 1.2, independently on the rel- 
evance of other evolutionary processes. 

To summarize, our study supports that major merg- 
ers have been the main drivers of the evolution of the 
massive end of the Red Sequence since z ~ 1.5, although 
other processes can also have contributed to it significantly 
at intermediate-to-low masses (especially, since z ~ 0.6). 
Our tests support observationally a late definitive buildup 
of the massive E-SO's through major mergers (mostly at 
0.6 < z < 1.2), in agreement with the expectations of hi- 
erarchical models of galaxy formation. 



9 CONCLUSIONS 

We study the buildup of the Red Sequence by analysing the 
structure, morphology, and star formation properties of a 
sample of red galaxies with stellar masses M,>5x 1O 1O M 
at 0.3 < z < 1.5. The novelty of this study is two-fold: 
first, our red galaxy sample includes galaxies both on the 
Red Sequence and at nearby locations on the Green Valley 
to trace transitory evolutionary stages towards it, and sec- 
ondly, we have simultaneously considered structural (regu- 
lar/irregular), morphological (disks/spheroids), and star for- 



mation enhancement (HSF/LSF) to define galaxy classes 
that can be directly associated with intermediate-to-late 
stages of major mergers, as well as with normal E-SO's and 
red regular disks. The redshift evolution of the fractions and 
number densities of each red galaxy type have been derived. 
Finally, these data are used to carry out a set of novel ob- 
servational tests defined on the basis of the expectations of 
hierarchical models, to test the hierarchical origin of massive 
E-SO's. 

Both the number densities and fractions of regular 
galaxies increase with cosmic time at 0.3 < z < 1.5, tracing 
the progressive buildup of massive normal early-type galax- 
ies during the last ~ 9 Gyr (mostly E-SO's, but also some 
red disks). The fractions of red regular spheroids (E-SO's) 
plus red disks increase by a factor of ~ 6 since z ~ 1.3 
down to 2 ~ 0.5, whereas the fraction of irregulars (major 
mergers) decreases by the same factor. However, the number 
density of red irregulars (major mergers) is constant down 
to z ~ 0.7, decreasing at lower redshifts. This means that 
the fraction of red irregulars decreases only relatively to that 
of regulars, pointing to their transitory nature, as already 
claimed in previous studies. We find that the number density 
distributions of regular spheroids (E-SO's) and of irregulars 
(major mergers) cross at z ~ 0.9 for M* > 5 x 1O 1O M0. 

Enhanced SFRs compared to the average SFR of the 
global galaxy population at each redshift at 0.3 < z < 1.5 
are hosted by most irregular disks in our red sample and 
by nearly half of the irregular spheroids (in agreement with 
the major merger-related nature of both types). At z < 0.7, 
enhanced SF can also be found in nearly all red regular disks. 
On the other hand, only ~ 25% of red regular spheroids (E- 
SO's) harbour enhanced star formation at all redshifts. 

The main result of this study is that we provide 
observational evidence pointing to the existence of a main 
evolutionary path among red galaxy types at 0.6 < z < 1.5, 
being: 

Irregular disk — > irregular spheroid — » regular 
spheroid. 

This track traces the conversion of blue disks into 
passive E-SO's through major mergers and dominates the 
buildup of the Red Sequence at z > 0.6, in excellent agree- 
ment with the expectations of the EM10 model. Our data 
are coherent with the prediction of this model about that 
the bulk of massive E-SO's with M* > IO^Mq have been 
definitively assembled through major mergers during the 
~ 2.2 Gyr peri o d ela psed 0.6 < z < 1.2, as also proposed 
bv lllbert et all < |2010h . 

Our results support observationally several expecta- 
tions of hierarchical theories of galaxy formation concerning 
the buildup of massive red galaxies: 

(i) Data is compatible with the fact that the massive red 
regular galaxies at low redshifts derive from the irregular 
ones populating the Red Sequence and its neighbourhood 
at earlier epochs up to z ~ 1.5. 

(ii) The progenitors of the bulk of present-day massive red 
regular galaxies seem to have been disks that have migrated 
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to the Red Sequence mostly through major mergers at 0.6 < 
2 < 1.2 (these mergers thus starting at z < 1.5). 

(iii) The formation of E-SO's that end up with M. > 
lO n M0 at z — through gas-rich major mergers seems 
to have frozen since z ~ 0.6. 
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